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ABSTRACT 

Context. The galaxy number counts has been traditionally used to test models of galaxy evolution. However, the origin 
of significant differences in the shape of number counts at different wavelengths is still unclear. By relating the most 
remarkable features in the number counts with the underlying galaxy population it is possible to introduce further 
constraints on galaxy evolution. 

Aims. We aim to investigate the causes of the different shape of the if-band number counts when compared to other 
bands, analyzing in detail the presence of a change in the slope around K ~ 17.5. 

Methods. We present a near-infrared imaging survey, conducted at the 3.5m telescope of the Calar Alto Spanish-German 
Astronomical Center (CAHA), covering two separated fields centered on the HFDN and the Groth field, with a total 
combined area of ~ 0.27deg'^ to a depth of if ~ 19 (3cr,Vega). By combining our data with public deep Jf-band images in 
the CDFS (GOODS/ISAAC) and high quality imaging in multiple bands, we extract if-selected catalogs characterized 
with highly reliable photometric redshift estimates. We derive redshift binned number counts, comparing the results in 
our three fields to sample the effects of cosmic variance. We derive luminosity functions from the observed K-hand in 
the redshift range [0.25-1.25], that are combined with data from the references in multiple bands and redshifts, to build 
up the K-hand number count distribution. 

Results. The overall shape of the number counts can be grouped into three regimes: the classic Euclidean slope regime 
{dlogN/dm 0.6) at bright magnitudes; a transition regime at intermediate magnitudes, dominated by M* galaxies 
at the redshift that maximizes the product 0*^^; and an a dominated regime at faint magnitudes, where the slope 
asymptotically approaches -0.4(q-|-1) controlled by post- A/* galaxies. The slope of the 7f-band number counts presents 
an averaged decrement of ~ 50% in the range 15.5 < K < 18.5 (d\ogN/dm ~ 0.6 — 0.30). The rate of change in the 
slope is highly sensitive to cosmic variance effects. The decreasing trend is the consequence of a prominent decrease of 
the characteristic density 4>*K.obs (~ 60% from z = 0.5 to z = 1.5) and an almost flat evolution of Af^ ^bs (1'^ compatible 
with = —22.89 ± 0.25 in the same redshift range). 
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1. Introduction 

With the a dvent of the large ph otometric surveys (such as 
COSMOS - IScoville et al]l200d UKIDDS - iLawrence enH 
120071 etc), the coverage of the extragalactic number 
counts (hereafter, NCs) on multiple bandpasses has greatly 
improved. Benefiting from the new generation of wide 
area cameras and dedicated telescope facilities, these sur- 
veys have provided large galaxy samples with a sig- 
nificant improvement in efficiency (|Gabasch et al.l l2008t 
iMcCracken et alll2003[ ). Even the traditionally more prob- 
lematic NIR surveys have made considerable progress in 
covering areas close to a square degree, up t o limiting mag- 
nitudes hardly reachable a few years ago (jSimpson et al.l 
[200llConselice et al.l[2008l) . 

The NIR galaxy counts usually have been consid- 
ered a particularly useful method to constrain galaxy 
evolution and cosmology, providing a simpler, less bi- 
ased, test of galaxy evolution models. However, although 
the overall shape of the NCs is well defined, published 



counts still exhibit a considerable scatter, making it dif- 
ficult to narrow down the evolution in a small mag- 
nitude interval. Additionally, the apparent simplicity of 
the NCs hides a mixture of evolving galaxy properties 
that complicates the interpretation of the observed fea- 
tures in terms of a single physical origin. Thus, it is 
not surprising that the explanation for the flatness of 
the if-band counts (relative to the optical NCs) or the 
presence of a sharp break in the slope around 16 < 
K < 1% arc still a matter of debate (Gardner et alj 
19931: [Cri stobal-Hornillos et al. 2003; Ou adri ct al.. .20071 
lovino et al.ll2005l : iTemporin et al.ll2008f) . 

The change in the slope of the NCs is a direct con- 
sequence of galactic evolution. Indeed, any feature in the 
shape of the NCs is closely related to the luminosity distri- 
bution of galaxies at a given epoch. By disentangling the 
relative contribution from the luminosity functions (here- 
after, LFs) at different redshift ranges, we will be able to 
identify the driving force behind the shape of the NCs. 
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Therefore, we aim to refine our understanding of the K- 
band NCs in the Ught of galaxy evolution by reconstructing 
the NCs in the K-hanA from rest-frame LFs in the mul- 
tiple bandpasses probed by the observed X-band at dif- 
ferent epochs. We benefit from the substantial amount of 
work done on derivin g accurate LFs at different redshifts 
in the opticia.1 bands (IWol f et al.' 2003:'Blanton ct al.' 200S 



Marchesini et al 



I Cole et al 



Chen et al.l 120031 : lllbert et al. 2005; Gabasch ct al. 



2001 



2007 ) and NIR bands dK ochane k et al 



iPozzetti et al.l 12 003: Fe ulner et al 



Cirasuolo et al.H2007t lArnouts et al. 200';^ to disen- 



tangle the multiple galactic populations at different redshift 
that assemble together to create the observed NCs, and to 
discern the true nature of the reported features on that 
distribution. 

The purpose of this paper is two-fold. First we present a 
NIR survey conceived to serve as reference for future spec- 
troscopic follow up. Second, we use these datasets and high 
quality photometric redshifts to derive NCs and LFs in the 
observed if-band, allowing us to analyze the underlying 
galaxy population responsible for the observed shape of the 
K-h&nd NCs. The paper is organized as follows: In fj^lwe 
describe the observations, data reduction and the compila- 
tion of complementary data sets. Sections [3] and |4] show the 
multi-wavelength characterization of the ii'-selected sam- 
ples and the procedure to derive photometric redshifts. In 
fj5] we investigate the reliability of these photometric red- 
shift estimates. Then, we use these if-selected samples in 
inland ^ presenting if-band NCs in three different fields 
and making use of the photometric redshifts to explore the 
redshift distribution and to probe the impact of cosmic vari- 
ance in the NCs. In ^we discuss the connection between 
NCs and LFs, and present our LF estimates along with pre- 
viously published results in multiple bands and redshifts. In 

we derive the expected NCs distribution from the cor- 
responding LFs. Section i}TOl summarizes the conclusions of 
the paper. 

Throughout this paper we use Vega magnitudes unless 
noted otherwise and adopt the current standard cosmology 
Hq = 70 km-is-iMpc"\ = 0.3 and f^A = 0.7. 
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Fig. 1. Sky coverage maps of the NIR datasets in the 
Groth and HDFN fields. The outlines represent. Top- 
The CAHA-fl' K'-band flanking fields at both sides of 
the WHT-INGRID JT-band coverage of the original Groth 
strip (ICristobal-Hornillos et aI]|2003D . and the mk J-band 
pointings of IVillar et al. I (I2008ir The J-band central im- 
age is shown for comparison with the HDFN if -band im- 
age covering the same area (~ 225arcmin^). Bottom - The 
CAHA-172k X-band coverage of the HDFN. The GOODS- 
HST footprints are also shown for comparison. 



2. NIR observations 

-ftT-band observations were carried out on two sepa- 
rate fields: one covering the flanking fields of the orig- 
inal Groth strip (Groth et al. 1994) at a(J2000.0) = 
IAHT^A?,', (5(J2000.0) = 52°28'41" (hereafter, Groth- 
FF), and the second centered at the Gre at Observatories 
Origi ns Deep Survey north field (HFDN; iGiavalisco et al.l 
[200l . a(J2000.0) = 12''36™49^ (5(J2000.0) = 62°12'58". 
Additionally, we made use of complementary infrared imag- 
ing which include s a full mapping of the Groth s trip in the 
J and Ks bands (jCristobal-Hornillos et al.]|2003l : hereafter, 
CH03) , and J-band obs ervations in the Groth and HDFN 
fields (jVillar et al.ll2008l : hereafter V08), consisting of three 
pointings covering at least ~ 60% of the K surveyed area 
in Groth, and a single pointing in HDFN completely over- 
lapping with the K exposure. Details of the observations 
are shown in Table [TJ Fig. [1] shows the layout of the two 
areas covered. 

The Groth-FF infrared observations in the if'-band 
were obtained using OMEGA-PRIME (hereafter, W), 
mounted on the prime focus of the 3.5m telescope at Calar 
Alto Spanish-German Astronomical Center (CAHA) dur- 



ing three runs in 2000 May 15-17, 2002 March 28-31 and 
2002 August 19-23. The K' filter was preferred to Ks for 
a better removal of thermal background. No significant dif- 
ference in {K' — Ks) galaxy photometry was found (see 
Section [^ . 

The il' camera was equipped a Ik x Ik HgCdTe Rockwell 
array that provides a field of view of 7' x 7' with 0.396"/pixel 
scale. In the first runs, a small misalignment on the sec- 
ondary mirror caused the PSF to be slightly asymmetric. 
This problem was enhanced in the 2002 run, causing the 
FWHM to be noticeably higher, and reducing the depth 
of the observations. The seeing conditions were generally 
average, with the FWHM ranging from 1.10" to 1.50" ex- 
cept for a few pointings with a higher seeing ('^1.9"). The 
background emission in the NIR is bright, non- uniform and 
highly variable. Therefore, in order to perform an accurate 
sky subtraction, sky dithered exposures are required. We 
used a hexagonal dithering pattern of 20" side, repeated 12 
times, shifting the center 3" each time to create the final 
exposure. After excluding the edges (having a lower cov- 
erage), each pointing in the Groth-FF covers 47 arcmin . 
The integration times were chosen to keep the count levels 
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within the linear regime of the detector. Typical exposure 
times were 31 x 2 s per dithered position (co-adds x indi- 
vidual exposure time), rejecting the first frame each time to 
avoid charge persistency problems. Total integration times 
were around 84 minutes per pixel. Note that the integra- 
tion time for the last run was increased to alleviate the 
alignment problem. 

The whole Groth-FF mosaic comprises 22 pointings of 
about 7' X 7', for a total of 822 arcmin^. 80% of the 
total area is above the average exposure time per frame 
(~ 80 min), with relatively shallow limiting magnitudes of 
if ~ 19.1 - 19.8 (3cr). The CH03 Groth strip observations 
covered ~180 arcmin^ in eleven 4' x 4' pointings, carried 
out with the INGRID instrument and iiT-band filter in the 
William Herschel Telescope, to a deeper limiting magnitude 
of if ~ 20.5 (3(7). The V08 central Groth strip observations 
covered 702 arcmin^ in three 15'xl5' pointings taken us- 
ing the OMEGA 2000 {n2k) instrument at 3.5m GAHA to 
a limiting depth of J ~ 22.5 (3a). The combination of the 
Groth NIR observations leads to a total area of 955 arcmin^ 
in the iC-band and at least 462 arcmin^ simultaneously ob- 
served in both the J and K bands. 

The HDFN infrared observations in the Ks-hand were 
performed with the CAHA 3.5m telescope using il.2k, which 
is an improvement over the previous fi' instrument. The 
detector is a 2kx2k HgCdTe HAWAII-2 Array, with a 
0.450"/pixel and a 15.4' x 15.4' field of view. Observations 
were performed over the course of one observing run in 2006 
May 15-17. The seeing conditions were slightly better than 
in Groth-FF, ranging from 0.9" to 1.1". The typical expo- 
sure times at each dithered position were 50 x 3 s. The total 
exposure time is 1.7 h. After excluding the noisy edges, we 
obtained a final image of 232 arcmin^ to a, K ^ 19.5 depth 
(3cr). The same area is covered in V08 J-band observations 
of HDFN to a limiting magnitude J ~ 21.9 (3f7). 

2.1. Data reduction 

The data were reduced in a standard way using a combi- 
nation of the UCM NIR reduction software (Cardiel et al.) 
and the IRAF^ software package XDIMSUM. The basic 
methods are outlined below. 

First, an average dark image with the same exposure 
time and number of co-adds is subtracted from the science 
images. Second, each science frame is flat-field corrected 
and the background emission is removed. Finally, all the 
frames are combined using an iterative method. 

The skyflat image employed in the second step is built 
from the science frames by combining all the images with a 
median filter to remove sources. The background emission 
image is created for each science image using the median 
combination of the same pixel in the 3 previous images 
and the 3 next images. The positions of several stars are 
used to determine the relative shifts between background- 
subtracted images. The images are aligned to a common 
reference using integer pixel shifts, to preserve the Poisson 
nature of the noise. Since the PSF is well sampled, we do 
not expect an increase in the average seeing. Object masks 
in the combined image are constructed using SExtractor 
(jBertin fc Arnoutill996[ ). In a second step (and following 
iterations), the construction of the skyflat and the back- 



^ http://iraf.noao.edu/ 



3 

ground estimation is repeated (and improved), this time 
masking out the sources detected in the previous step. 

Several additional reduction procedures are carried out 
to improve the quality of the final images. We create a 
mask of bad pixels in each image. An initial cosmetic defect 
mask is created using the dark images. We then inspect each 
background subtracted image individually; images with se- 
vere artifacts, significantly higher seeing or very poor trans- 
parency are discarded, while others with localized artifacts 
(e.g. satellite trails) are masked using a custom procedure. 
Additional bad pixels in each image are identified using 
a cosmic ray detection procedure or are removed with a 
sigma-clipping algorithm during image combination. 

2.2. Photometric calibration 

Photometric calibration was performed comparing aperture 
photometry from Tw o Micron AU Sky Survey (2MASS) 
(jSkrutskie et al.ll2006l ) bright stars in the final mosaic af- 
ter rejecting those with poor quality flags in the 2MASS 
catalog. Despite the different K-hand filters, the integrated 
Kg fluxes are essentially preserved, as the color transfor- 
mations (Eq[T|) are typically smaller than a few 0.01 mag, 
due to the smooth behavior of NIR SEDs. The n'-K' and 
2MAS S-ji'g magnitudes are rela ted by the following equa- 
tions (jWainscoat fc Cowielll993 ): 

Ks = K + Q.005{J-K) (1) 
K = K' - (0.22 ± 0.03)(i7 - K). 

We find that the same equations apply for Vl2k with no 
significant dispersion. 

The number of stars employed for each field calibration 
ranges from 7 to 15, leading to a zero-point rms between 
0.04 and 0.08 in Groth-FF, and - 0.05 in HDFN, which 
should only be considered as a lower bound to cumulative 
uncertainties introduced by the color transformations. To 
check the quality of the Groth-FF photometric calibrations 
we compared them to CH03 INGRID K-hand observations 
inside an ^80 arcmin^ ove rlapping region. Using standard 
stars (jPersson et al.l [r998l) . they found their photometric 
uncertainty to be less than 0.03-0.05 mag. The median off- 
set of the comparison is less than 0.07 mag for bright ob- 
jects, so we do not attempt to readjust the zero-points to 
their calibration. 

Note that the shift in effective wavelength between the 
K' and Kg filters is small, and much less than the filter 
widths. In the discussion that follows we do not distinguish 
between the different K filter sets. 

2.3. Catalog completeness 

We estimated the completeness of our catalogs by simulat- 
ing the detection and photometry of fake sources. The fake 
sources were created by extracting a bright source from the 
image, scaling it to the desired flux level, and injecting it at 
random locations in the central well-exposed regions of the 
images. We then attempt to detect these objects by run- 
ning SExtractor under the same parameters as in the origi- 
nal frame. Figure [2] shows the resulting completeness curves 
for point sources as a function of magnitude for some of the 
pointings. Because the simulated sources are not required 
to fall on empty regions, the confusion to real sources may 
slightly lower our detection efficiency. Nonetheless, given 
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(J2000.0) 
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(s) 


(arcscc) 


(mag) 


(mag) 


[arcmtn ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


grothll 


14:14:43.515 


+52:03:14.76 


4800 


1.15 


19.1 


19.3 


47.2 


grothl2'* 


14:15:28.298 


+51:57:35.62 


5820 


1.53 


18.7 


19.3 


48.7 


groth21 


14:15:11.712 


+52:08:13.31 


4680 


1.20 


18.8 


19.2 


47.1 


grot/i22 + 


14:15:56.989 


+52:02:48.05 


7500 


1.50 


18.8 


19.3 


45.4 


grothSl 


14:15:39.456 


+52:13:11.31 


3720 


1.50 


18.7 


19.2 


47.4 


groth32 


14:16:21.756 


+52:07:41.58 


2520 


1.30 


18.7 


19.1 


46.6 


groth41 


14:16:06.669 


+52:18:38.36 


3240 


1.16 


18.3 


18.5 


47.2 


groth42 


14:16:49.262 


+52:12:49.23 


5040 


1.31 


19.1 


19.4 


48.6 


grothSl 


14:16:35.898 


+52:23:32.60 


5040 


1.00 


19.1 


19.5 


49.6 


groth52 


14:17:17.680 


+52:18:02.82 


5040 


1.14 


19.2 


19.5 


48.5 


groth61 


14:17:07.477 


+52:28:30.29 


4680 


1.30 


19.2 


19.7 


48.9 


grotfi62 + 


14:17:50.812 


+52:22:59.62 


6180 


1.50 


18.7 


19.3 


47.0 


groth63 


14:17:04.981 


+52:23:23.18 


5040 


1.14 


18.8 


19.1 


49.1 


groth64 


14:17:17.456 


+52:30:05.59 


5040 


1.07 


19.2 


19.5 


46.8 


groth71 


14:17:35.943 


+52:33:22.61 


5040 


1.42 


18.9 


19.6 


50.9 


groth72 


14:18:21.179 


+52:28:19.52 


4920 


1.80 


19.0 


19.8 


47.4 


grothSl 


14:18:03.943 


+52:38:41.66 


5040 


1.30 


19.2 


19.8 


47.4 


groth82 


14:18:49.704 


+52:33:22.23 


2940 


1.90 


18.4 


19.3 


48.2 


grothgi 


14:18:33.323 


+52:43:46.63 


4560 


1.50 


19.2 


19.7 


48.7 


groth92 


14:19:17.913 


+52:38:09.21 


4920 


1.60 


18.9 


19.5 


50.7 


grothlOl 


14:19:02.220 


+52:48:38.00 


4800 


1.90 


18.5 


19.7 


49.4 


grothW2'* 


14:19:45.792 


+52:43:11.56 


6300 


1.50 


18.9 


19.5 


49.6 


gooodsn 


12:36:49 


+62:12:58 


6300 


1.1 


19.2 


19.5 


232.2 



Table 1. Units of right ascension are hours, minutes and seconds, and units of decHnation are degrees, arcminutes and 
arcseconds.(f) Pointing observed during secondary mirror misaUgnment problem. Col. (7) is the 5a limiting magnitude 
measured inside a 1" radius circular aperture. 



the uncrowded nature of our K-hand images, we expect 
the confusion to affect our completeness values by less than 
5% at the - 80% level. 

For a more realistic depth estimation, we have simulated 
extended sources as well as point like sources. Extended 
sources are considered when the effective radius is 20% 
greater than the typical radius of the stars on that frame. 
As expected, the efficiency for extended sources lowers the 
estimation from point-like sources by 0.2 to 0.3 magni- 
tudes. A precise determination of the completeness correc- 
tion would require measuring the detection efficiency over 
a wider range of effective radius. Nevertheless, for the pur- 
pose of our investigation, we can establish a conservative 
magnitude threshold of K=18.5 in the HDFN and Groth- 
FF without affecting our conclusions. Hence, excluding the 
shallowest fields (groth32, groth41, groth82, grothlOl) from 
the final catalog, we obtain a galaxy sample of reasonably 
homogeneous depth in both fields, avoiding the use of com- 
pleteness corrections. The estimated 80% completeness and 
limiting magnitudes for each pointing are given in Table [1] 
For the Groth strip data, CH03 carried out simulations find- 
ing insignificant completeness corrections below K = 19. 
We will adopt that value as the limiting magnitude. 

2.4. Complementary data 

In addition to our NIR survey, high-quality imaging and 
photometry are publicly available for both fields. For the 
Groth field, we made use of some of the panchromatic data 
sets that have been acquired as a part of the All- wavelength 
Extended Groth Strip International Survey (AEGIS, see 
iDavis et al.ll2007l for a detailed data description), including 
ground based ugriz deep imaging from the Canada-France 
Hawaii Telescope Large Survey (CFHTLS-D03-, Gwyn et 
al. in preparation), observed with MegaCam at the 4 m 
CFTH; BRI wide field (0.70° x0.47°) observations using 
the CFHT12K mosaic camera; Spitzer mid-IR data cov- 
ering the wavelength range 3.6/iTO to 8.0/im (GTO pro- 



Band 






Source 


(1) 


(2) 


(3) 


(4) 


IRAC-3.6.. 


3.561 


21.6 


Spitzer GTO 


IRAC-4.5.. 


4.510 


21.8 


Spitzer GTO 


IRAC-5.8.. 


5.689 


21.8 


Spitzer GTO 


IRAC-8.0.. 


7.958 


21.7 


Spitzer GTO 


U 


0.358 


25.2 


Subaru deep imaging'^ 


B 


0.442 


25.2 


Subaru deep imaging"' 


V 


0.546 


24.9 


Subaru deep imaging" 


R 


0.652 


24.4 


Subaru deep imaging" 


I 


0.795 


23.9 


Subaru deep imaging" 


z 


0.909 


23.6 


Subaru deep imaging" 


b 


0.430 


25.7 


GOODS'' 


V 


0.592 


24.9 


GOODS'' 


i 


0.770 


24.3 


GOODS'' 


z 


0.906 


23.9 


GOODS'' 


HKs .. 


2.127 


21.3 


QUIRC deep imaging" 


b,spectra.. 


0.430 


24.4 


TKRS" 



Table 2. Main characteristics of the datasets in HDFN. 
Col. (1): Observing band. Col. (2): Effective wavelength 
of the filter calculated by convolving the Vega spectrum 
(Colina & Bohlin 1994) with the transmission curve of the 
filter -|-detector. Col. (3): Limiting AB magnitude. Col. (4): 
Source from where the data were obtained. ° Publicly 
available ultra-deep optical and NIR data from Capak et 
al.(2004). The Great Observatories Origins Deep Survey 
(GOODS; Giavahsco et al. 2004a). " Team Keck Treasury 
Redshift Survey (TKRS; Wirth et al. 2004) and Cowie et 
al. (2004). 



gram) and a deep R-b and image from Subaru-SuprimeCam 
(jMivazaki et aLlboOTl ). 

In the HDFN, we used ultra-deep optical and NIR data 
spanning from th e U to the HK-ha,nd {UBVRIzHKs, 
ICapak et ani2004f) . together with HS T-ACS bviz imagin g 
published by the GOODS Team (Gi avalisco erani2004[ ). 
The main characteristics of each data set are given in 
Tableland Table E] 
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Band 






Source 


(1) 


(2) 


(3) 


(4) 


Kp 


2.114 


20.3 


CAHA-Oprirne 


R 


0.652 


- 


Subaru deep imaging 


u 


0.381 


27 


CFHTLS 


g 


0.486 


28.3 


CFHTLS 


r 


0.626 


27.5 


CFHTLS 


i 


0.769 


27 


CFHTLS 


z 


0.887 


26.4 


CFHTLS 


B 


0.437 


24.5 


CFHT-12k 


R 


0.660 


24.2 


CFHT-12k 


I 


0.813 


23.5 


CFHT-12k 


J 


1.255 


22 


WHT-INGRID 


K 


2.170 


21.8 


WHT-INGRID 


IRAC-3.6 


3.561 


21.6 


Spitzer GTO 


IRAC-4.5 


4.510 


21.8 


Spitzer GTO 


IRAC-5.8 


5.690 


21.8 


Spitzer GTO 


IRAC-8.0 


7.957 


21.7 


Spitzer GTO 


nuv 


0.232 


24.5 


GALEX GTO 


fuv 


0.154 


24.5 


GALEX GTO 


U 


0.361 


24.8 


INT-WFC 


B 


0.436 


25.5 


INT-WFC 


R,redshift.... 


0.660 


24.2 


DEEP2 



Table 3. Main characteristics of the datasets in Groth. 
Col. (1): Observing band. Col. (2): Effective wavelength 
of the filter calculated by convolving the Vega spectrum 
(Colina & Bohlin 1994) with the transmission curve of the 
filter+detector. Col. (3): Limiting AB magnitude. Col. (4): 
Source from where the data were obtained. 



At the time of writing, the CFHT12k fully reduced im- 
ages were not publicly available. Instead, raw images and 
calibration files were downloaded from the CFHT-CADC 




Fig. 2. The K-h&wd best fit to the point source com- 
pleteness curves. Point-like sources were inserted at ran- 
dom locations in the images. The completeness is defined 
as the fraction of recovered sources in each image. The 
short dashed line shows the completeness curve for the 
Groth32 pointing, the shallowest frame employed comput- 
ing the number counts. The vertical long dashed line at 
K = 18.5 depicts the magnitude threshold to which our 
catalog is ~ 100% complete. 



archive and the Elixir website^, and were reduced using 
the IRAF MSCRED package. Precise astrometry and pho- 
tometry calibrations were carried out by comparison to the 
CFHTi2k pubhc catalogs, and the CFHTLS and Subaru 
imaging in overlapping regions. 

To complement the imaging data we have compiled a 
set of spectroscopic redshifts obtained by several surveys 
on these field s. For the HDFN w e have 1699 spectros copic 
redshifts fro mlWirth et all (|2004f ). ICowie et all (|2004f ). and 
iReddv et al.l(|2006D . Most of these sources are below z 1 
and have a high reliability flag 80%). In the Groth 
field we have ^15000 redshifts from the DEEP2 collabora- 
tion ()Davis.2006 ). Only a small fraction of those are found 
within our surveyed area. 

3. Source characterization 

3.1. Multi-wavelength photometry 

Multicolor photometry in all available bands was ob- 
tained employing a similar method to the one described in 
Perez-Gonzalez et al. (2008, Appendix A). Briefly, the K- 
band source catalog, obtained with SExtractor, was cross- 
correlated in a 2" radius to each one of the UV, optical and 
NIR catalogs. Sources detected in fewer than 3 bands were 
rejected. Given the relatively shallow limiting magnitude of 
the if -band images, fully covered by deep optical and IRAC 
images, this is a safe method to reject spurious sources. 
Then, aperture matched photometry was performed using 
elliptical Kron-like apertures. The aperture size and orien- 
tation is determined by the K-hand image and translated 
to all the other optical and NIR images except for the IRAC 
bands. 

For our typical seeing values the aperture is large 
enough to enclose the PSF in all these bands. For the 
IRAC bands, whe re the resolution is slightly worse (^^2"; 
iFazio et al.l l2004f ). the flux was measured in small cir- 
cular apertures (typically 3") and corre cted using stel- 
lar PSF growth cu rves (similarly to iHuang et al.1 12003 : 
iBarmbv et al.ll2008h . In addition, we have also apphed a 
deblending algorithm for the sources with multiple if -band 
counterparts associated with a single IRAC source and sep- 
arated by more than 1" (the approximate limit of the as- 
trometric resolution). The if -band positions of the sources 
were used to re-align the position of the aperture in the 
IRAC band. Then, the IRAC PSF is convolved with the K- 
band PSF and the flux is measured again inside a 0.9" aper- 
ture, applying the corresponding aperture correction. For 
source separations larger than 1", t he flux contamination 
with this method is lower than 10% (|Perez-Gonzalez et al.l 
12001 . 

For GALEX and HST data, we took the SExtractor 
MAG_BEST magnitude of the closest source. Consequently, 
we do not use these data in the photometric redshift deter- 
mination. For the rest of the bands, the procedure allows 
us to obtain accurate colors, since we measure a similar 
fraction of the flux from an object in each filter. 

Uncertainties in the measured flux were derived taking 
into account the background noise, photon statistics, read- 
out noise and scatter in the WCS. The standard method 
of determining the background noise from pixel-to-pixel 
variations of adjacent sky pixels often underestimates the 

^ http:/ /www. cfht.hawaii.edu/Instruments/Ehxir/ 
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real value due to correlat ed signal effects introduced dur- 
ing t he data reduction (|Labbe et al.l l2003UGawiser et all 
l2006t) . To obtain a more accurate determination of the 
background noise we follow ed the method described in 
Perez-Gonzalez et alj (j2008l ). which is similar to that of 



Labbe et al. I (|2003D . For each source, the sky flux was mea- 



sured in randomly distributed apertures with the same size 
as the photometric aperture. Then, we estimated the back- 
ground noise from the width of the flux histogram, approx- 
imated by a Gaussian distribution. The comparison to the 
sky background measured with SExtractor shows that the 
latter tends to underestimate (10%-15%) the noise even in 
the K-hand images, where integer pixel shifts were used 
during the frame combination. 

3.2. Photometric redshift 

We calculate photometric redshifts using the methods de- 
scribed in detail in Perez Gonzalez et al. (2005,2008). 
Briefly, a reference set of galaxy templates is fitted to the 
observed spectral energy distributions (SEDs), taking into 
account the flux uncertainties. The galaxy templates are 
composed of a sample of galaxy spectra from HDFN and 
CDFS with highly reliable redshift determinations, includ- 
ing some z > 1.5 galaxies, and well-covered SEDs (with 
more than 10 different photometric data points). The set is 
fltted with models of single, and composed stellar popula- 
tions (1-POP and 2-POP) and models of dust emission. The 
comparison between data and models is done using a maxi- 
mum likehood estimator that takes into account the uncer- 
tainties in each data point. The resulting 1-POP templates 
are characterized by four parameters: The star formation 
time scale t, age t (assuming exponentially declining star 
formation history), metallicity Z , and extinction A(V). For 
the 2-POP we have twice the same family of parameters 
to characterize both the young (instantaneous burst) and 
the old populations. The final set comprises 3624 models 
(1666 1-POP-hdust, 1958 2-POP+dust). The photometric 
redshift estimation and uncertainty is derived from the min- 
imization of the probability distribution obtained from 
the fitting of the observed data to a grid of redshifted mod- 
els (using S z — 0.01). In additio n, we include AGN tem- 
plates from iPolletta et aT] (|2007D to provide a better fit of 
the very few AGN dominated SEDs. 

Furthermore, it is possible to use the best fitting SED to 
estimate rest-frame fiuxes, luminosities and also observed 
fluxes in bands that might be missing due to shallow imag- 
ing or differences in the covered area. The use of these syn- 
thetic magnitudes allows a consistent calculation of colors 
for all sources in the sample. 



3.3. Star-galaxy separation 

Stars are separated from galaxies using a combination of 
morphological and color criteria. Following Perez- Gonzalez 
et al. (2008), all secured objects (with detection in more 
than 3 photometric bands) are classified as stars if they 
satisfy at least 3 of 9 morphological and color criteria. The 
first one is the all-band weighted SExtractor STELLARITY 
parameter. When available, the STELLARITY estimations 
from HST images are given a higher weight in the com- 
parison. The other criteria are based on nIR and IRAC 
color-color criteria and comparisons between magnitudes 



deriv ed from isophota l and circul ar apertures extracted 
from Eisenhardt et al.1 (|2004f) and iRowan-Robinson et al.l 
(|2005|) (see Perez-Gonzalez et al. 2005, 2008 for a de- 
tailed description) . We also consider the BzK criteria by 
iDaddi et al.1 (|2004[ ). to isolate the stellar locus: {z — K)ab < 
0.3 X {B — z)ab — 0.5. In Fig. [3] we compare the efficiency 
of the latter to the other color- morphology criteria. Sources 
classified as galaxies are required to satisfy fewer than three 
different criteria. We find that almost all the stars satisfy 
the BzK criteria with only a few objects below our star 
threshold, probably being galaxies scattered into the stel- 
lar boundary due to photometric err ors as already reported 
by other authors (|Kong et al.ll2006^ . 

We also classify as stars a very small fraction of objects 
whose best fit to a spectral template is the Vega SED, and 
that does not have a reliable spectroscopic redshift. These 
objects typically satisfy at least 2 of the other criteria, and 
after visual inspection, it is most likely that they are stars. 





B-z 



Fig. 3. Efficiency of the BzK color-color diagram to iden- 
tify stars, compared to the combination of the other 9 cri- 
teria applied to isolate stars. The dashed lines depict the 
regions defined by the BzK criteria. The bottom line iso- 
lates the stellar locus, while the other two select sources at 
z>1.4. Black stars show the sources that satisfy three or 
more stellarity criteria. 



4. Composite sample 

Our final K-selected sample is comprised of 1313 sources in 
the 223 arcmin^ of HFDN, 4660 sources in the 648 arcmin^ 
of the Groth-FF, and 1957 sources in the 161 arcmin^ of 
Groth strip. Out of these, less than 15% are identified as 
stars in each field. We merged both Groth samples into a 
single catalog with 5948 sources. To avoid repeated sources, 
the two samples were cross-correlated in a 1" radius remov- 
ing duplicate sources with lower S/N. 

Additionally, to complement the low depth sam- 
ples, we have created a K-band selected sample in 
the Chandra Deep Field South (CDFS, a(J2000.0) = 
03''32™28^ (5( J2000.0) = -27°48'30"). The reference im- 
age was drawn from the latest data release of the 
GOODS/ISAAC observations (Retzlaff et al. in prepara- 
tion). This final release adds two new frames in the K-band, 



Barro et al.: On the nature of the extragalactic number counts in the Tf-band 



7 



HDFN 
-CDFS 
- Groth 




-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 

AZ/(1+Z) 




0.2 0.4 0.6 0. 



photo 



Fig. 4. Comparison of the scattering: in the photometric „. _ t-.i , . • i i -n i- , -i ■■ r .^ ^ ^ 

J , .p, ,. , / \ /fi ? \ / ii Fig. 5. Photometric redshiit distributions oi the A-band 

redshift estimates (z.pec - Zphot)/(l + ^phot) for the three ^^^^^^^^ ^^^.^^ ^^^^ ^^^^ ^^^^^ ^^^^^ 

nelds. ° 



increasing the total surveyed area to 172 arcmin^. The com- 
bined mosaic covers the region with a variety of exposure 
times, with an average depth of if ~ 22.7 and a maxi- 
mum of i^T 24.2. Consequently, we can safely consider 
the whole sample to be complete bel ow K = 20 (see also 
the area-depth estimations from Graz ian et an[2006 for the 
vl.5 data). The CDFS has also been the target of extensive 
multi- wavelength observations. We m ake use of all the pub- 
lic and proprietary data compiled in lPerez-Gonzalez et al.1 
(|2005f ) to create the multicolor sample, following the proce- 
dures described in section ^ After excluding a small por- 
tion of the mosaic with very low S/N regions, the sample 
contains 6810 galaxies in an area of 168 arcmin^. 

5. Photometric redshift distribution 

We have checked the quality of our photometric redshifts 
with all the available spectroscopic data compiled for the 
three fields {Fig.\^. For HDFN there are 287 sources (36% 
of the sample) with highly reliable spectroscopy. The av- 
erage redshift difference {Sz = Zspec — Zphoto) is 0.001, 
70% of the sample present values of (Jz/{l -I- z) < 0.05 
(where cr^ is the absolute value of Sz), and 91% have 
(7z/{l + z) < 0.1. In the CDFS 232 sources (53% of the 
sample) have spectroscopic redshift with a high quality flag. 
The mean value of Sz is 0.014. 72% of the sample present 
values of a:,/{l + z) < 0.05 and 91% have cr^/(l + z) < 0.1. 
In the Groth-FF we find 846 sources with high quality flag 
(35%). The average Sz = 0.036. 59% of the sample present 
values of cr^/(l + z) < 0.05 and 85% have (7^/(1 + z) < 0.1. 
Finally, there are 334 spectroscopically measured sources 
in the Groth strip with an average Sz — 0.010. The 68% 
and 89% have values of crz/{l + z) below 0.05 and 0.10 
respectively. 

Figure [5] shows the photometric redshift distribution for 
the samples in Groth, HDFN and CDFS to a limiting mag- 
nitude of K = 18.5. The distributions have been derived 
taking into account the typical photometric redshift error. 
Although the shape of the distributions are typical for a flux 
limited sample, field-to-field differences reveal the effects of 
large scale structure. The peak of the redshift distribution 
depends on the luminosity function of the galaxies at dif- 



ferent redshifts. At K < 18.5 the majority of the sources 
are found at z < 0.7 (see next section), leading to an expo- 
nentially decreasing tail at higher redshifts. The prominent 
density peak in HDFN at z ~ 0.5 is in good agreement 
with the redshift distribution found by the Team Keck spec- 
troscopic survey IjWirth et al.ll2004 ). The less pronounced 
feature in CDFS at approximately the same redshift also 
coincides with a spectroscopic ally confirmed overdensity at 
z = 0.7 (|Vanzella et al.ll2006D . if we take into account the 
broadening caused by the photometric redshift uncertain- 
ties and the decreasing selection function at high redshift. 
Finally, there appears to be a slight underdensity at z ~ 0.4 
and a small peak at z ~ 0.2 in the Groth field. This alterna- 
tion of peaks and voids between fields highlights the impact 
of field-to-field variance in the small area surveys, especially 
at z < 1. 

6. K-band number counts 

Figure [S] shows the if-band NCs for the three fields, as 
well as the NCs from CH03, for comparison. The NCs in 
0.5 magnitude bins, up to ^=18. 5 for Groth and HDFN, 
and up to K—20.0 for CDFS are summarized in Table [H 
Within these limits, detection efficiencies for point sources 
are well above 90% (see Fig. [21 and completeness estimates 
in ij2.3l and Q. Hence, we do not apply completeness cor- 
rections to the NCs. Our counts do not need corrections for 
spurious sources either; spurious sources typically appear 
at brightness levels where the efficiency drop s significantly 
below 100% fCH03! Eliche-Moral et al.ll2006l i. and the frac- 
tion of spurious sources in our catalogs must be negligible 
given that each source is required to be detected in at least 
three bands f ^3T|) . NCs are tabulated in Table H The er- 
ror calculation assumes Poisson statistics for low numbers 
(jGehrelsl Il986f ) added in quadrature to the standard de- 
viation derived from bootstrapping the source magnitudes 
convolved with a Gaussian error kernel. 

In the Groth field, our counts are in good agreement 
with those of CH03. Differences are Icr compatible in the 
range 16.0<K< 18.5. The overall agreement provides an ex- 
ternal check on the quality of both sets of if-band number 
counts in the Groth field. As we move to brighter magni- 
tudes, our counts are systematically above those of CH03 
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Fig. 7. Left: Averaged if-band number counts from Groth, HDFN and CDFS compared with a compilation of results 
taken from various sources. Right: Slope for the differential number counts in our three fields, along with data from the 
literature. The horizontal dashed lines indicates the value of the slope at both sides of if = 17.5 reported in CH03. 
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Table 4. Differential number counts in 0.5 magnitude bins 
for the X-selected samples in HDFN, CDFS and the Groth 
field, uncorrected for completeness. Only the NCs up to the 
completeness limit of each field are shown. 
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Fig. 6. The raw number counts in our three fields (HDFN, 
CDFS, Groth) and counts from CH03 obtained in an small 
portion of our total area in the Groth field. Error bars are 
derived from a combination of Poisson errors and boot- 
strapping. No completeness corrections have been made. 

and HDFN. Such an offset is likely due to the smaller area 
mapped in these two surveys (~5 times smaller in the case 
of CH03). Nevertheless, NCs below if < 16 probe very low 
redshift populations (see next section) that would require 
larger areas to be studied accurately. 

The NCs in the Groth and HDFN fields at K > 18.5 are 
shown in Fig.[6]for comparison to the deeper CDFS counts 
and the completeness corrected counts of CH03. The fact 
that the shallower NCs do not fall abruptly until K ~19 
supports our efficiency estimates for the Groth and HDFN 
images. In addition, it is noteworthy that the scatter in the 
brighter counts {K < 17) is on average greater than the 
typical uncertainties for the NCs in any of the fields. The 
most plausible cause of this discrepancy is cosmic variance. 
As we showed in the previous section (Fig. [51), there might 



be significant overdensities at some redshifts that lead to 
these differences. Additionally, the effect is accentuated by 
the fact that bright counts come mainly from low redshift 
sources, which are poorly sampled in small area surveys, 
like HDFN and CDFS. We will further discuss field-to-field 
variance effects in the next section. 

The left panel of Fig. [7] shows the averaged NCs 
of the 3 fields along with counts from the literature. 
We have compiled measurements from large area sur- 
veys to mitigate the effects of cosmic variance on the 
shape of the NCs. That is th e case of the 1 . 47 de g^ of 
the DEEP2/Palomar survey (IConselice et al.ll2 008). the 
0.70 dcg2 of the UKIDDS -UDS dSjm pson et alfeooel the 
0.17 deg2 of the Daddi-F (|Kong et al. 2006]) and the com- 
bined -0.11 deg^ of the MUSYC survey (jQuadri et all 
2007). Additiona lly, we compared our re sults with the 
deeper NCs of iForster Schreiber et al.l (|2006D in the 
MS1054-03 galaxy cluster to illustrate the asymptotic be- 
havior of faint counts and also the impact on the NCs of a 
very prominent overdensity such as a cluster at z=0.83. It 
can be seen in Fig. [7] that the density peak dominates the 
bright counts. This is comparable, to a lesser extent, to the 
effect at if < 17 of the low redshift peak in CDFS. 
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The logarithmic derivative of the differential NCs (i.e., 
the slope of the NCs) is shown in the right panel of Fig. [T] 
The horizontal lines at 0.25 and 0.6 indicate the value of 
the slope on both sides oi K = 17.5 reported in CH03. 
The thick black line shows the evolution of the slope for 
the combined NCs of our three fields. We confirm that 
the logarithmic slope decreases rather sharply from 0.6 at 
K = 16.25 to ~0.3 at K =18.0, a similar resuh to CH03. 
However, it appears that the evolution of the slope is best 
described by a continous decreasing trend, rather than by 
fixed values on both sides of ~ 17.5; i.e., the slope of K- 
band NCs decreases monotonically over the 16. 0< K <19.0 
range. The logarithmic slopes derived from counts by other 
authors are bound by the same upper and lower limits as 
in our data, except at K < 16.0. However, they show a 
more gentle variation with magnitude, as well as a greater 
dispersion. Such variations among different surveys may 
be driven by two processes. First, completeness and re- 
liability of the photometry are treated with varying de- 
grees of rigor by different authors, which contributes to 
the scatter in the NCs (see CH03). Second is cosmic vari- 
ance. Substructure in the redshift distribution can cause 
significant fiuctuations in the characteristic density, and 
hence slope variations, that might lead either to a sharp 
break or to a smoother evolution in the slope of NCs. 
These density peaks are clearly recognizable in the right 
panel of Fig. [7] at A' = 16.5 and K 18 when compar- 
ing the NCs of HDFN and CDFS. To strengthen this idea 
we h ave created 100 Millennium simulation mock catalogs 
from iKitzbichler fc Wtiit^ (|2007D selected over a 0.25 deg^ 
area, similar to that of our Groth field, to sample the extent 
of fie ld-to-field variance effects. A lthough the if-band NCs 
from IKitzbichler fc Whitj (j2007[ ) do not accurately repro- 
duce the observed distribution a.t K > 17 (mainly because 
the z > 1 galaxy population is overestimated compared to 
the observations), their results on clustering analysis are 
consistent with observations. Thus, the simulation is suit- 
able to mimic the effects of cosmic variance. As can be seen 
in Fig. [U the broadening effect in the slope of the bright 
counts is similar to ours, whereas the confidence interval for 
the faint counts is significantly narrower, probably because 
at that point our NCs are dominated by photometric errors 
and depth effects, and not by large scale structure (LSS). 

7. Redshift distribution of K-band number counts 

The change in the slope of NCs is an indirect effect of galac- 
tic evolution. Indeed, any feature in the shape of NCs is 
closely related to the luminosity distribution of galaxies at 
a given epoch. By disentangling the relative contribution 
to the NCs of the LFs at different redshift ranges, we will 
be able to identify the main effect responsible for shape of 
the NCs. 

Fig. [9] depicts the if-band NCs in the CDFS, the deepest 
sample. We have used the photometric redshift estimation 
to split the total number counts into redshifts bins. The 
binned NCs tend to resemble the shape of a (Schechter) 
LF. The evolution of the comoving distance with redshift 
causes an effective shift in the successive redshift bins to- 
wards fainter magnitudes, and the posterior accumulation 
of the higher redshift NCs. In addition, despite the likely 
greater effects of cosmic variance on a single field, it seems 
clear that for K > 18 the total number counts tend to 
become a mixture of very different redshifts ranges, with 




Fig. 8. Slope for the differential number counts in the 100 
samples of the Millennium simulation d rawn from the mock 
catalogs of IKitzbichler fc White! (|2007D . The thick dashed 
line shows the averaged NCs combining all the samples. 
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Fig. 9. NCs in the if-band for the CDFS field alone (thick 
black line) and NCs in several photo-redshift ranges. At 
bright magnitudes {K^17) the contribution to the NCs is 
limited to only a few low redshift bins, whereas at if ~18 
the total number counts becomes a mixture of several low 
and high redshift bins {z > 1.5) 



a growing contribution of high redshift galaxies (z >1.5) 
absent at brighter magnitudes, whe re 2-3 low redshift bin s 
account for 90% of the total NCs (jConsehce et al.ll2008D . 
Furthermore, it can be seen that for i^T < 16 the counts 
sample the bright part of the LFs at z <1, causing values 
to be highly sensitive to both volume effects and LSS. 

The four panels of Fig. (TUj show the relative contribu- 
tion to NCs in HDFN, CDFS and the Groth field divided 
in four redshift bins: [0.25-0.50], [0.50-0. 75], [0.75-1] and [1- 
1.25]. Below z < 0.25 the averaged contribution is small 
enough (< 10%) to be neglected. The next three bins are 
responsible for most of the total counts up to K ^ 18. 
However, the exact proportion might differ by up to a factor 
of 2 due to LSS. The alternate peaks and valleys between 
HDFN/CDFS and the Groth fields between z = OA ~ 0.9 
(see Fig. [5]) cause the most prominent differences around 
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K = 17.0 — 17.5 in the first panel of Fig. (TUl similarly to 
the underdensity in CDFS at z ~0.8 which leads to the 
significantly lower NCs around K ~ 17.5 — 18 in the third 
panel. 

8. Number counts & luminosity functions 

In the previous section, we showed evidence for a mild de- 
creasing trend in the slope of the NCs. In this section, 
we will try to find the origin of such a decreasing trend 
by studying the LP at several redshift intervals (benefit- 
ing from the estimation of photometric redshifts for the 
entire sample). The calculation of LFs must take into ac- 
count that the observed K-hand probes progressively bluer 
rest-frame bands with increasing redshift. By z ^ 0.75, the 
/C-band central wavelength shifts to 1.25/im( J-band), and 
for z > 1.5 it begins to probe the optical bands. Therefore, 
in order to derive NCs from LFs it is necessary to select 
different rest-frame bands at each redshift bin. 

In iJSll, we present the function al relation betwe en the 
NCs and LFs, parametrized using a lSchechterl (|1976( ) func- 
tion, that we will use to derive the if-band NCs. In 5Hl2, we 
derive LFs from the observed if-band, in the redshift range 
[0.25-1.25], using our three galaxies samples. In SJ513 we 
combine our LFs with LFs from the literature to explore the 
general picture of luminosity evolution in the optical and 
NIR. Finally, in iJ5l4, we summarize the multi- wavelength 
LFs probed by the observed i^-band at different redshifts, 
that gives rise to the isT-band NCs distribution. 

8.1. Number counts from Schechter functions 

The distribution of observed galaxy counts is a consequence 
of the LFs and the cosmological framework. Hence, any 
feature in the NC distribution could be explained in terms 
of the evolution of the LFs assuming a cosmological context. 
Here we will adopt a ACDM framework, and assume that 
the galaxy LF can be described by means of a Schechter 
function, 

0(A/) = 0.41n(10)(/)*10°-4(*^*-*^)("+i) exp (lO°-4(*^*-^'^)) (2) 

where M* is the characteristic absolute magnitude, a the 
faint-end slope, and (j)* the normalization of the luminos- 
ity function. Although a Schechter function might not be 
appropriate to fully describe the faint end population of 
the LF (Blanton et al. 2005), this is not the case for the 
bright NCs, which will be mostly dominated by M* galax- 
ies (see next section). The NCs represent the distribution of 
galaxies per apparent magnitude and sky area. Assuming a 
parametrization for the LF, the functional form of the NCs 
becomes 

f^f dVc 

'N{mx„,\M,(j),a\ xp ) = / (j)(m,z,\M,(j),a] aq )—^dz{3) 
(1+^) J^. (1+^) ail 

where Zi and Zf are the lower and upper limit of the red- 
shift bin, dVc/dQ is the differential comoving volume, and 
m\M,(f),a\ \„ represents the Schechter parameters of a 

LF in a band with effective wavelength ttt-t- This indi- 
cates the explicit dependence of the NCs in the ii'-band on 
bluer LFs with increasing redshifts. Note that even if we 
assume that the Schechter parameters are constant in the 



redshift bin {zi,Zf], the expression [3] still depends on LFs 
at different wavelengths. When Az is small, we can approx- 
imate the NCs from single LFs in a photometric band at 
■pq^, where z is the mean value of the redshift bin (z^, Zf]. 

Additionally, Eq. [3] can be interpreted as a sum of LFs 
at different redshift bins, weighted by the corresponding 
comoving volumes. Hence, the shape of the NCs would be 
the result of the LF parameters at a given epoch and their 
evolution with redshifts modulated by the volume element. 
Furthermore, the slope of the total NCs is the sum of the 
slope from each redshift bin weighted by the normalized 
NCs. 



/ dN \_ 1 djV _^Ni/ 1 dNiX 
^^\dmdn) N dm ^NKNidm) ' 

where Ni are the NCs in the redshift bin i, and N is the sum 
of all the redshift bins. The slope of the NCs presents some 
interesting properties. First, it is independent of the abso- 
lute value of (j)* , but it does depend on the relative change 
of this parameter (i.e., it depends only on the evolution of 
4>*), which can be normalized arbitrarily. Second, in addi- 
tion to the low redshift Euclidean limit {d\ogN/dm =0.6), 
it can be shown that, at very faint magnitudes, the slope 
approximates asymptotically to — 0.4(q; -|- 1). 

\dmd\lJ m» ^I{a,M{m)) ' 

where /(a;,M(m)) accounts for all other terms except 
]^g-o.4m^ including the integral over redshift and the ex- 
ponential term in the Schechter equation, and M(m) in- 
dicates the relation between absolute and apparent mag- 
nitudes through the luminosity distance. The /(a, M{m)) 
in the numerator differs from the denominator only in the 
derivative of the exponential term, which also outputs the 
2Q-o.4m ^g]-jjj^ indicated explicitly. At faint magnitudes the 
first term of the Eq. [S] tends to zero, dominated by the 
faster decrement of 10"° '*'", yielding the asymptotic limit 
— 0.4(a -|- 1). More intuitively, it can be understood as a 
regime where the NCs are mostly dominated by the faint 
end of a single LF. Hence, the slope of the NCs necessarily 
becomes the slope of this LF. For any given variation with 
redshift of the Schechter parameters, there will always be a 
maximum value of the product (j)* ■ The LF at that red- 
shift will be favored among the others (suppressed smaller 
values of the product) and, at sufficiently faint magnitudes, 
it will be the main contributor to the NCs. 

A very interesting result that emerges from the analysis 
of the NCs in the full magnitude range is that the evolu- 
tion of the slope can be summarized in three main regimes, 
namely: 

— (1) The Euclidean regime. The classical low red- 
shift approximation, that yields the well-known result 
dlogN/dm = 0.6. The NCs in this regime are mostly 
populated by M< M* at very low redshift (z< 0.2). 

— (2) The transition regime. The slope departs from the 
Euclidean limit into a weighted sum of the slopes of LFs 
at low-mid redshifts. The dominant contribution would 
come from the LF at the redshift that maximizes the 
product 4'*jf (the factors controlling the normaliza- 
tion of Eq.[3li. In this regime, the slope decreases rapidly 



Barro et al.: On the nature of the extragalactic number counts in the Tf-band 



11 





Fig. 10. Relative contribution to the K-hand number counts in rcdshift bins for the samples in the HFDN, CDFS and 
Groth field. Field-to-field differences can reach the 40% for the most prominent LSS features. 



around the apparent magnitude of M* at the dominant 
LF (i.e. at the knee of the LF). In the absence of a sig- 
nificant evolution in </)* , the maximum would take place 
at th e peak of the volume element at z~2 (jBershadvl 
120031) . 

(3) The "a regime". In this regime, the majority of the 
LFs dominating the NCs contribute with galaxies at the 
faint end. Therefore, the slope will be a combination of 
individual slopes approaching the minimum value (^- 
0.4(a+l)). As a consequence of the previous regime, 
this phase should be controlled by the same LF as in 
the previous phase. However, other LFs, not strongly 
suppressed by the (t>*^f factor, and having a signifi- 
cantly larger a, might dominate this phase. 



We will further discuss the evolution of the slope in terms 
of this schema in iMl. 



8.2. LFs from the observed K-band 

In order to derive NCs in the ii'-band using Eq. [3] we need 
to feed the equation with LFs in different bands at differ- 
ent redshifts. More precisely, we are interested in the LFs 
from the observed K-hand. Thus, we have obtained LFs in 
several redshift bins at decreasing rest-frame wavelengths 
centered at 2.16/im/(l -I- z), where z is the mean redshift 
of the bin (z^, zj]. 



8.2.1. Methodology 

To esti mate the observ ed LF we have applied the 1/Knax 
method ()Schmidtlll968D. We associ ate the Poisson errors to 
each magnitude bin ( Gehrelsl[l986l ) added in quadrature to 
photometric errors. Since we have not applied k-corrections, 
Am translates directly into AM. Additionally, in order to 
quantify the uncertainties due to the photometric redshift 
errors, we have performed Monte Carlo simulations at each 
redshift bin, assuming an uncertainty characterized by the 
width of a Gaussian distribution derived fitting the his- 
togram of 0-2/(1 + z) for each redshift bin. 

Additionally, we measured the LF from the STY 
(jSandage et al.ll979l ) estimator, which is a parametric max- 
imum likelihood method. The STY is not as sensitive to 
large-scale fluctuations as the 1/Vmax method, but requires 
a priori fixing the functional form of the LF. Nevertheless, 
the assumption that (j){M) is well described by a Schechter 
function suits our purpose of explaining the NCs in terms 
of LFs, as we have described in the previous section. 

To check the effect of photometric redshift errors in the 
STY method, we perform the same Monte Carlo simula- 
tions as for the 1/Vmax- Since the probability of galax- 
ies existing at redshift z grows rapidly until z ^ 2.5, a 
larger fraction of intrinsically faint sources are shifted to 
the bright end than viceversa. This can potentially bias the 
determination of the LF introducing large systematic effects 
when Az is l arge a nd the distance m odulus grows rapidly 
()Chen et all l2003l : iMarchesini et all [2007il . Nevertheless, 
the combination of the typical values of az ~ 0.05 — 0.07 
with the relatively shallow limiting magnitude of our LFs 



12 



Barro et al.: On the nature of the extragalactic number counts in the Tf-band 




Fig. 11. LFs derived from the observed /C-band in different redshift slices. The LF with the Vmax method are shown 
as open squares with error bars (including Poisson, photometric and photometric redshift uncertainties). The continuous 
line depicts the best fit from the STY method. The inset shows the best fitting (M,a) values (star) with 1,2 and 3(7 
confidence contours, and also the mean value (cross) with error bars derived from the photometric redshift simulations. 
The vertical line shows the photometric threshold for the shallowest fields, -ft^voga = 18.5. The additional vertical line in 
the first panel shows the lower threshold at bright magnitudes. 



(less than ~ 1 mag deeper than M*) cause only a moderate 
offset in the Schechter parameters (AM < 0.2,Aa < 0.1). 
It can be seen from the inset of Fig. [11] that the mean and 
standard deviation resulting from the simulations (star) lie 
always within the la-2a confidence level derived from the 
best fitting result of the STY method. 

8.2.2. Observed luminosity functions 

The LFs, (j){M)dAI, of the combined samples were com- 
puted in four redshift bins of width Az ~ 0.5 using the 
methods described above. 



redshift bin 


Mit,o5.(AB) 


C^K.obs 


rK.otsiW-'Mpc-') 


0.25 - 0.50 
0.50 - 0.75 
0.75 - 1.00 
1.00 - 1.25 


-22.95(fixed) 
-23.0lg;li 
-23.08S!:i| 
-22.96g:i 


-i.oogii 
-i.02g:i? 

-0.93g:i 


n.a.iU 
•3-41.70 
Q ^0.22 
•J-^O.IS 
Q nO.n 
J.Zq 24 
q0.31 
^••JO.26 



Table 5. Best fitting Schechter parameters for the STY 
LFs from our combined samples. The la errors were derived 
from Monte Carlo simulations 



Table [5] summarizes the best fitting values of the 
Schechter function for each redshift. Fig. [TT] shows the esti- 



mates from the T4riax method (open squares) with Icr error 
bars (from Poisson statistics and photometric redshift sim- 
ulations), and the LF from the maximum likelihood method 
at different redshifts. The inset shows the 1, 2 and Scr con- 
fidence levels for the a and M* parameters (open star) , to- 
gether with the mean value derived from the photometric 
redshift simulations (cross). The vertical line indicates the 
photometric threshold of the shallowest sample (i4r=18.5). 
For fainter luminosities, the values of the LFs are derived 
exclusively from the CDFS sample. Note that the com- 
pleteness limit of our samples translates into luminosities 
very close to or even fainter than M* for the higher red- 
shift bins. This means that galaxies fainter than M* will 
not contribute significantly to the NCs up to if < 18.5. 
Nevertheless, the deeper coverage of the CDFS allow us to 
properly fit the LFs that will be employed to derive the NC 
function (see section H8A\\ . 

In order to properly account for cosmic variance, the 
values of for the Vinax method derived exclusively from 
the CDFS data (those above the photometric threshold; 
vertical line of Fig. fTTj) have been corrected by a scale factor. 
This scale factor has been derived from the median value 
of all fields in the magnitude bins in common, weighted by 
the area covered by each field. Also, only the magnitude 
bins with more than 100 galaxies in the CDFS sample have 
been considered in the calculation. 
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Finally, for the 0.25 < z < 0.50 bin, the bright extremity 
of the LF is poorly constrained even in the total combined 
area of the three fields. Nevertheless, since the faint end 
is properly sampled, we decided to fit the LF fixing the 
value of M* and setting a threshold on the bright {K < 16) 
magnitudes. The value of M* was taken from the literature. 
At z ~ 0.38 the rest- frame wavelength probed by the K- 
band is close to the iJ-band. However, no references for 
iJ-band LFs at that redshift were available then, and we 
decided to use M* from the res t -frarn e if-band LF at z ~ 
0.4 published bv lArnouts et al.l (|2007[ ).Then. we corrected 
that value by applying an H-K color term. The mean value 
of that color at z ~ 0.4 derived from our data is {H — 
K)(AB) = —0.21 ± 0.11, very similar to the local value 
(~-0.20; iJarrett et al1l2003[ ). 

8.3. LFs evolution in rest-frame bands 

Analyzing the evolution of the LFs directly derived from 
the observed K-hsmd requires exploring LFs in shorter rest- 
frame wavelengths at different redshifts. For this purpose 
we compiled LF estimates from the literature in optical and 
NIR bands. We have co r npiled the l ocal LFs in the J a nd K - 
bands from ICole et al.l (|200l[ ) and iKochanek eFall (|200l[ ) 
respectively, whereas for the iJ-band we estimated a value 
of M* by applying a color correction to the K-hand value, 
following the same procedure as in the previous section. 
The local optical LFs in the i,z and R bands were drawn 
from lBlanton et all (|2003f ). 

At higher r edshift , the iiT-band LFs were drawn from 
lArnouts et all (|2007( ). Their LFs for the SWIRE- VVDS- 
CFHTLS survey predicts an evolution in the Schechter pa- 
rameters which is consistent with previous results from 
ICirasuolo et all |2007) in the ukidds -UDS. For the J-band 
we considered the LF estii nates fromlPozzetti et"aLl ()20Q3l ) 
(K20 survey, 52 arcmin^V H^ulner et aLl (|2003h (MUNICS 
survey, 0.17 deg^), and iDahlen et al] (pOOl (GOODS- 
CDFS, 130 arcmin^). Their results are compared in Figure 
15 of the latter, showing that, despite the apparent inconsis- 
tency in the Schechter parameters, there is good agreement 
in the data points. The different observational constraints 
of each survey lead to different best fitting values for a dis- 
tribution with a very similar sha pe. We have ta ken the LF 
in th e J-band at z ~ 0.48 from lFeulner et aTl ( (2003. ) and 
from IDahlen et all (|2005) at z ~ 0.9. The lar ger area of 
the first is better to constraint the bright end at low red- 
shift, while the depth of the latter is more suitable at higher 
redshifts. 

We find the same apparent inconsistencie s due to the 
a — M * degeneracy in the LF s of optical bands. Illbert et al.l 
(|2005D (VVDS, 0.6 deg^) and lGabasch et al.l ([20061) (FORS 
Deep Field Survey, FDF, 35 armin^) derived very differ- 
ent values of the Schechter parameters that nevertheless 
produce 1 — 2a compatible LFs when comparing both es- 
timates under a common limiting magnitude and using 
Monte Carlo simulations (jCabasch et al.l I2006D . For our 
purposes we are more interested in the high-redshift op- 
tical LFs (z > 1 .5). T hus, we preferred the estimates from 
iGabasch et al.l (|2006D that, despite the smaller area, arc 
able to probe the LF 3-4 magnitudes deeper. In addition 
to the LFs at z = 1.60 and z = 2.26, we derived some 
values for M* and (p* at other redshifts using the evolu- 
tionary parametrization also given in Gabasch et al. (2006; 
as (1 + z)'^). Finally, the LF in the R-band at z 2.3 was 



drawn from iMarchesini et al.1 (|2007f ). The authors derive 
high redshift LFs in optical bands from a combination of 
very deep nIR infrared observations (including 51arcmin^ 
to if ~ 23). Table ini summarizes the Schechter parameters 
of the LFs in the different bands and redshifts. 

The two panels of Fig. [T^j show the evolution with 
redshift of M* and (j)* in several NIR and optical bands. 
Vertical dotted lines indicate the effective rest-frame wave- 
length of each band. The open circles show the values of 
M* and (p* in each band at different redshifts increasing 
from top to bottom. For clarity, the local values have been 
plotted with small black dots. The dashed lines connect 
the values of the Schechter parameters for the different 
bands at the same redshift (continuous grey line at z=0). 
Finally, the large black dots show our estimates of the LFs 
at z '--^ 0.4,0.6,0.9,1.1. These estimates do not match ex- 
actly any standard photometric band. However, the value 
at z ~ 0.4 is almost centered in the _ff-band, and the LFs 
at z ~ 0.6 and z ~ 0.9 should be consistent with the esti- 
mates in the J-band. Comparing our values with the results 
for these bands of Dahlen et al. and Feulner et al., we find 
a very a good agreement in the characteristic luminosity, 
whereas the values of (j)* are poorly consistent. A possible 
explanation for this discrepancy is that the values of Dahlen 
et al. at z ~ 0.9 might be affected by an underdensity peak 
between z = 0.7 — 1.1, while our results present a slight 
overdensity around that redshift (see Fig. [5]). Nevertheless, 
the comparison to the values of Feulner et al. is 2a compat- 
i ble, a nd is also consistent with the results of lPozzetti et al.l 
(|2003l ). Wider areas are necessary to reduce the field-to- field 
variations and obtain more reliable estimates. 

It can be seen from the left panel of Fig. [T^j that the 
multi-band values of M* at z = resemble the shape of a 
galactic SED, with the absolute magnitude peaking around 
the H band, close to the 1.6/xm stellar bump. At higher 
redshifts, the shape is not preserved due to the different 
luminosity evolution in each band; i.e., in the K-hand, a 
substantial brightening of ~ 1 mag between the local uni- 
yerse and z ^ 2 has been reported by several authors 
("C aputi et all 120061 : ICirasuolo et~an l2007t ISaracco et all 
200^. In the J-band, Sar acco et al. find AM* ^ -0.7 to 
z ~ 3, and the results from'Marche sini et al.l (|2007l ) indicate 
a similar or slightly higher brightening (AM* ~ 1.7) in the 
optical LFs up to the same epoch. However, it is not surpris- 
ing that the brightening is faster in the optical bands than 
in the NIR since they are more affecte d by the light com- 
ing from younger ste llar populations (jDahlen et al.l [20051 : 
Illbert et al.l[2005t Imif et al...200a'l . 

In the right panel of Fig.[T2[we can see that the charac- 
teristic density, </>*, follows a decreasing trend in every band. 
The average estimates for optical LFs indicate a decrease of 
^ 50% in the number density from the local value to z ~ 1 
(Gabasch et al. 2006, Ilbcrt et al. 2005), whereas the evo- 
lution in the NIR LFs seems to be sli ghtly weaker (^ 30%) . 
Als o, the results in the Jf- band from lArnouts et al] ([20071 ) 
and ICirasuolo et al] ([2007D suggest a mild decreasing trend 
from z = 0.4 — 1.25. Finally, the LFs of IMarchesini et al] 
(2007) in optical bands confirm the decreasing trend at high 
redshift finding a decrement of the 70 — 80% in 0* from 
the local value to z~ 3. 
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Fig. 12. Evolution with redshift of the Schechter parameters M*(left panel), 0* (right panel) in the RizJK rest-frame 
bands. For each photometric band (RizJK, vertical dotted lines) the empty circles depict the redshift evolution of the 
parameter at 2; 0.4, 0.9, 1.3 and 1.9, drawn from different authors (see Table [5]). The long dashed lines connect the 
values of the parameter at the same redshift in different bands. The black stars shows the value of the parameter at 
z — 1.70, z — 2.26 and z = 2.25 in the i, z and R bands respectively. The black dots indicate the LF parameters at 
z = 0.38, 0.62, 0.88, 1.12 derived in this work. 




Fig. 13. Multi- wavelength evolution of the Schechter parameters Af^^f,^ (left panel), 4>*k obs (right panel) as probed by 
the observed if-band. For each photometric band (RizJK, vertical dotted lines) the empty circles depict the redshift 
evolution of the parameter at z ^ 0.4, 0.9, 1.3 and 1.9, drawn from different authors (see Table[5]). The black stars show 
the value of the parameter a,t z — 1.70, z — 2.26 and z — 2.25 in the i, z and R bands respectively. The grey vertical line 
at z = 1.25 is the estimated redshift from which the contribution of the sources to the NCs up to K — 18.5 is less than 
10%. The black dots indicate the LF parameters at z = 0.38,0.62,0.88, 1.12 derived in this work. The thick black line 
shows the approximated evolution of the Schechter parameters as probed by the observed iiT-band at growing redshifts. 



8.4. LF evolution from the observed K-band 

The two panels of Fig. [13] shows the multi- wavelength evo- 
lution of the Schechter parameters as probed by the ob- 
served if -band (thick black line); i.e, LFs at progressively 
shorter rest- frame wavelengths for growing redshifts. This 
is essentially the same plot as in Fig. [T^] except that now 
the X-axis indicates the redshift where the if-band probes 
the different rest-frame bands. The vertical grey line cen- 
tered at z = 1.25 indicates the approximate redshift limit 
where the contribution to the K-hand NCs {K < 19) be- 
comes smaller than 10%. Althought this means that the 
impact of the optical LFs on the bright counts is negligible. 



they illustrate the process of construction of NCs from LFs, 
allowing us to estimate the NCs to fainter magnitudes. 



Note that in the optical bands, we have no references 
for the LFs at the precise redshift. Therefore, together with 
the estimated value of the z-band LF at z = 1.42 (open 
circle), derived from the Gabash et al. parametrization (see 
previous section), we show the measured value at z = 1.7 
(black star) , which presents an M* significantly fainter than 
their estimates at lower redshift. The same applies for the 
LF in the i-band at z = 2.26 (instead of z = 1.89) and the 
R-band at z — 2.25 (instead of 2 = 2.36). 
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Source 




Limit 

[mag] 


Area 


Band 


z 


Af * — 5 log hro 
[AB] 


r 


a 


Kochanck ct a] 


(2001) 


K< 11.25 


< 7000 


K 


0.02 


-22.36 ± 0.05 


3.97 ± 0.34 


-1.09 ± 0.06 


Ainouts ct al. 


(2007) 


F(3.6Aim) > 9tJ.Jy 


0.85 


K 


0.50 


-22.83 ± 0.30 


3.47 ± 1.4 


-1.1 ± 0.2 


Arnouts ct al. 


(2007) 


F(3.6Aim) > g/ijjf 


0.85 


K 


0.90 


-23.12 ± 0.08 


4.00 ± 0.3 


-1.11 ± 0.06 


Arnouts ct al. 


(2007) 


F(3.6tim) > 9tJ.Jy 


0.85 


K 


1.35 


-23.13 ± 0.17 


2.60 ± 0.5 


-1.1 ± 0.2 


Ainouts ct al. 


(2007) 


F(3.6/jm) > 9tJ.Jy 


0.85 


K 


1.75 


-23.28 ± 0.14 


1.62 ± 0.3 


-1.1 ± 0.2 


Fculncr ct al. 


(2003) 


K< 17.50 


0.18 


J 


0.30. .0.60 


-22.93 ± 0.24 


2.60 ± 0.80 


-1.00 fixed 


Dahlcn ct al. 


(2005) 




0.036 


J 


0.75. .1.00 


-23.09±»i* 


1 Q7+0.60 




Gabasch at al. 


(2006) 


I(AB)< 26.8 


0.01 


z 


0.4 


-22.78 ± 0.41 


2.7 ± 1.2 


-1.33(fixed) 


Gabasch at al. 


(2006) 


I(AB)< 26.8 


0.01 


z 


0.9 


-22.93 ± 0.44 


2.2 ± 0.7 


-1.33(fixed) 


Gabasch at al. 


(2006) 


I(AB)< 26.8 


0.01 


z 


1.32 


-23.05 ±0.48 


1.9 ±0.5 


-1.33(fixed) 


Gabasch at al. 


(2006) 


I(AB)< 26.8 


0.01 


z 


1.31. .1.91 


-22.71 ± 0.18 


2.0 ±0.3 


-1.33(fixed) 


Gabasch at al. 


(2006) 


I(AB)< 26.8 


0.01 


z 


1.89 


-23.14 ± 0.51 


1.70 ± 0.4 


-1.33(fixed) 


Gabasch at al. 


(2006) 


I(AB)< 26.8 


0.01 




0.4 


-22.58 ± 0.31 


2.7± 1.3 


-1.33(fixed) 


Gabasch at al. 


(2006) 


I(AB)< 26.8 


0.01 




0.9 


-22.69 ± 0.34 


2.2 ± 0.8 


-1.33(fixed) 


Gabasch at al. 


(2006) 


I(AB)< 26.8 


0.01 




1.32 


-23.78 ± 0.37 


1.9 ± 0.5 


-1.33(fixed) 


Gabasch at al. 


(2006) 


I(AB)< 26.8 


0.01 




1.89 


-22.85 ± 0.40 


1.7 ± 0.4 


-1.33(fixed) 


Gabasch at al. 


(2006) 


I(AB)< 26.8 


0.01 




1.91. .2. 61 


-22.93 ± 0.14 


1.9 ± 0.2 


-1.33(fixed) 


Marchcsini ct al. (2007) 


K < 23(max) 


0.65 


R 


2.0. .2. 5 


-22.67±!J;^o 


1 n7±*^-^'^ 


n-|+0.21 



Table 6. Summary of the derived parameters of the LF parameters obtained by various authors by fitting them with 
Schechter functions, in several rest-frame bands, in different redshift bins, and from different samples of field galaxies. 




O) 0.6 
E 



NCs encompass 
confidence interval (1o) 

- K-band 

- 4i=3.4; z>0.5 

- o[=-1.50; z>0.5 

- LFcte 




Fig. 14. Left: if-band NCs derived from the LFs (thick black line) compared to data from other authors (the legend as 
in the left panel of Fig. [7]). The shaded region indicates the la confidence interval derived from the simulations on the 
uncertainties of the Schechter parameters. The continuous lines with different symbols show the redshift binned number 
counts at z = 0.32,0.62,0.88,1.12 derived from the LFs in the present work. The da shed lines with different symbols 
sho w the redshift binned nu mber counts at higher redshifts from lGabasch et"an (j2006D (in the i and z rest- frame bands) 
and iMarchesini et al.l ()2007[ ) (in the R rest-frame band). These NCs have a negligible effect on the on the total NCs at 
i^<18. The arrows depict the best-fit value of M^^^^ for the LFs at ^ = 0.32,0.62,0.88,1.12. Note that to derive the 
total NCs we have also used LFs at higher redshift from the literature. Right: Slope of the NCs derived from the LFs 
(thick black line). The dark grey shaded region represents the la confidence interval derived from the simulations on 
the uncertainties of the Schechter parameters. The light grey shaded region shows the region encompassing the slopes 
derived from the references. The dashed line shows the predicted slope derived from a single LF constant through all 
redshifts. The circled line shows the predicted slope fixing a = —1.50 for the LFs at z > 0.5, The crossed line shows the 
predicted slope fixing (f)}^ = 3.4 x lO^^h^QMpc^^ for the LFs at z = 0.5. 



8.5. Summary 



In this section we presented the tools and measurements 
required to recover the NCs in terms of LFs at different 
bands and redshifts. As outlined at the beginning of this 
section, our motivation here is to determine the family of 
LFs at different redshifts that constitutes the NCs in the 
if-band, a result which is easily summarized in the thick 
black lines of Fig. [T31 



As it can be seen from the left panel of that figure, the 
characteristic luminosity ^j,^ shows an almost flat evo- 
lution beyond z — 0.4. In fact, it is la compatible with a 
constant evolution: = -22.89 ± 0.25. On the other 

hand, (I)k obs shows a progressive decline of the ~ 60% 
from the local value in the rest-frame if-band to the high 
redshift value in the R-band, with a significant decrease 
around z > 1. A number of recent results seem to confirm 
the existence of this decline, suggesting that it is mostly 
driven by a steeper decrement in the number density of 
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quiescent galaxies ()Zucca et alj 120061 : lArnouts et al.l l2007t 
iFaber et al.ll2007l ). This trend is usually interpreted as an 
indicator of the beginning of an epoch of major build-up for 
these populations. Nevertheless, a deeper analysis of the 
evolution with redshift of the color bimodality is beyond 
the scope of this paper. 

The evolution described by the {K, obs) quantities has 
no direct physical meaning regarding galaxy evolution. 
However, it is crucial to describing how galaxy evolution 
is encrypted in the i^-band NCs, and provides an inter- 
pretation for the shape of the NCs in terms of meaningful 
quantities, as we will show in the next section. 

9. iC-band NCs from evolving LFs 

Using the multi- wavelength evolution of q^^, (p'^ pre- 
sented in the previous section (thick black line of Fig. [T3|) 
and their corresponding values for the a parameter (see 
Table [5]), we have derived the NCs in the K-hand and 
its slope up to /iT < 19 (these results are summarized in 
Table [7]). Additionally, we have calculated the 68% confi- 
dence interval for both using Monte Carlo simulations on 
the Schechter function parameters. The random values of 
M* and a for each LP were taken from the error ellipses 
calculated with the STY method. When the ellipse param- 
eters were not available, such as for the optical bands, we 
used a Gaussian distribution with a median equal to the 
given value and a equal to the la deviation. 

In the left panel of Fig. [T3] we compare the if-band 
NCs derived from the LFs (solid black line) with the obser- 
vational results shown in the left panel of Fig. [7] (the leg- 
end is the same). The remarkably good agreement proves 
the validity of the method applied to relate both quanti- 
ties, and allows us to interpret the variation of the slope 
in terms of LFs. Therefore, we can consistently argue that 
the main causes of the shape of the iiT-band NCs are: the 
almost flat evolution with redshif of ^^j, and the signif- 
icant decline of (j)*^ ^^j, with redshift. The former causes the 
redshift binned NCs to become progressively concentrated 
at fainter magnitudes due to the variation with redshift of 
the distance modulus (higher at low redshift but decreasing 
rapidly at higher redshifts). However, the dominant effect is 
the decrease of the characteristic density, which opposes the 
growing comoving volume. The relative growth in the co- 
moving volume per redshift bin (Az = 0.25) progressively 
decreases from a factor of 2 at z=0.32 to 1.3 at z=1.12 
(as we approach the peak at z~2). Simultaneously, (px obs 
presents a decreasing trend that matchs the growth of the 
volume element at z=1.12, freezing the contribution of the 
LFs at higher redshifts, that otherwise would continue to 
increase to the peak of the volume element. 

In terms of luminosity, it can be seen in the left panel of 
Fig. [TUthat the dechne of ^bs leads to a smooth transi- 
tion between the intrinsically faint galaxies at z < 0.5 and 
the ~ obs galaxies at 0.5 < z < 1 around K ^ 17. b 
(the second and third arrows of the panel represent ^j,^ 
at those redshifts). Hence, since the slope of the LFs (and 
therefore of the NCs) decreases rapidly after the knee, the 
transition leads to more flattening compared to the slope 
at if = 16, where the ~ M* galaxies at z < 0.5 clearly 
dominate the NCs. 

To help us clarify the relevance of Af^^j,^ and <j)*i( obs 
for the modeling of the slope, in the right panel Fig. [HI we 
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Table 7. Theoretical if -band NCs and NCs slope in 0.5 
magnitude bins derived from the LF of section [S] 



show the slope of the if -band NCs resulting from assuming 
the local if -band LF through all redshift bins (constant LF; 
dashed line). Despite the poor resemblance to the slope de- 
rived from the observed LFs (i.e. , the slope of NCs in the left 
panel; solid black line), the decreasing trend is indicative 
that a gradual change in the slope is a natural consequence 
of the shape of the LFs modeled by the volume element. 
However, a precise evolution in the LFs is required to re- 
produce the change in the slope in the appropriate place. 

To strengthen this idea we have also explored the sce- 
nario of no density evolution, fixing the value of obs 
z > 0.5 (crossed line). In the absence of a decreasing char- 
acteristic density, the weigth of the M> M* galaxies up to 
z~2 would lead to a slope systematically larger than our 
la prediction at if > 17.5, and clearly out of the envelope 
of the observed NCs at if > 18.5. 

Finally, we have tested the role of a in the shape of the 
NCs. As can be seen from both panels of Fig.[T31 the rather 
flat values of a in the low redshift LFs accentuates the de- 
creasing trend of the slope. However, these values might be 
slightly underestimated due to the relatively shallow depth 
of the samples. The line with open circles in the right panel 
of Fig. [T4l shows the predicted slope with a = —1.50 fixed 
in the LFs at z > 0.5. As expected, the steeper faint end 
does not significantly affect the NCs around if ~ 17.5, 
which are dominated by M* galaxies. However, it does pre- 
dict a higher slope at fainter magnitudes that falls outside 
the region defined by the observations. This test poses a re- 
striction against high values of a > —1.5 in the low redshift 
(z=0.5-l) NIR LFs. 

9.1. Summary and predictions for fainter magnitudes 

Within the three regime schema for the evolution of the 
slope proposed in section H8.1\ it can be seen that the 
slope of the if -band NCs leaves the Euclidean limit around 
if ^ 15.5. It reachs the (M* dominated) transition regime 
faster (K~17.5) than in a no evolution scenario (K>19) 
due to a significant decrease in (j)'^ that freezes the con- 
tribution from bright z >1 galaxies. The migration into 
the a regime is harder to determine. The relatively low 
(a = —0.9, —1) faint end slope at z^l might cause LFs with 
larger values of a to dominate the faint counts. These LFs 
would not be severely suppressed if the decreasing trend 
in (f)'^ obs stabilizes at z > 1.5, given that the volume ele- 
ment evolves smoothly between z=1.5-4. Therefore, for any 
a <-l the high-z LFs could eventually dominate the NCs 
at sufficiently faint magnitudes. 
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For the set of LFs that we have compiled, Gabash et 
al. propose a fixed a =-1.33 for the LFs at z=1.3-2.6 in 
the i,z bands, whereas Marchesini et al. give a — —1.07 
in the R-band at z—2.35. Interestingly, several authors re- 
port signi ficantly larger values o f a for the optical LFs 
at high- z (|Marchesini et al.ll2007l a =-1.40 in the B-band 
at z=3; iReddv et all l2008l a =-1.80 at 1700l at z=3). 
Based on these numbers, we predict a slowly decreasing 
slope from d\ogN / dm ~ 0.3 — 0.2 in the magnitude range 
K =19.5-22, due to the dominant contribution from the 
z~2 bin (q!=-1.33). However, the slope in this range is 
not purely in the asymptotic a regime (that would lead 
to d\ogN / dm ~ 0.13). The resulting value is a weighted 
sum of slopes which includes a contribution from z > 2 
(around K ~ 20) M* galaxies. This contribution from 
larger slopes lowers the rate of approach to the asymptotic 
limit. Furthermore, if the slope of the high-z LFs grows 
rapidly to large values (a ~ —1.8) we could measure an 
increase in the slope around K >23. 



9.2. NCs in other NIR bands 

The procedure to derive NCs from LFs described in section 
[5] can be used to reproduce the NCs in any other band using 
the appropriate LFs and redshifts. Nevertheless, based on 
the multiwavelength LFs probed by the observed if -band 
(summarized in Fig[T3] for the if-band) it is possible to 
predict the approximate shape of the NCs in the closest 
NIR bands; e.g., IRAC-3.6 or H (1.65Atm). 

Firstly, if the NCs in different bands are referred to 
AB magnitudes, the results can be compared in a more 
meaningful context, avoiding the offset introduced by the 
AB-Vega transformation. In the AB-system, the magni- 
tude range where the NCs fall in the Euclidean regime 
{dlog N / dm ~0.6) is determined by the local value of M* 
and (f>* in that band. Therefore, similar values of the local 
LFs will result in similar NCs at bright magnitudes. 

Furthermore, if the multiwavelength evolution of M* 
and (j)* probed by the different bands exhibits a similar evo- 
lution, the slope of the NCs should change around the same 
magnitude range. That is the case of the H and K bands. 
At z=0, we have derived the H-band LF directly from the 
if-band applying a color term. At higher redshifts, the mul- 
tiwavelength LFs probed by the observed iJ-band ( J-band 
at z=0.32 and z-band at z=0.83) presents similar values 
to the LFs probed by the iiT-band at the same redshifts 
(see Fig. [T3]). As a consequence, the NCs in the i?-band 
also pres ent a flattening in th e slope around m[AB]'^19- 
Metcalfe et al.l [2 0061 and references therein 



In 



Magdis et al. 12008) have shown that the NCs 



19.5 (see 
addition, 

in the IRAC-3.6 band exhibit a significant flattening at 
slightly fainter magnitudes (m[3.6]^20). This suggests that 
IRAC-3.6 and K probe a similar evolution in the LFs. If we 
estimate the [3.6] local LF following the same approxima- 
tion as for the i?-band, the shift in the slope flattening is 
consistent with ap plying the typical color K-[3.6]'~-0.6 mag 
()Dale et al.ll2007[ i to the K-hand local LF. 

A detailed study of the change in the slope would de- 
pend on the multi-band LFs probed by the given band. 
Nevertheless, this simple calculation provides a rough esti- 
mate of the magnitude where the change in the slope begins 
to be appreciable. 



10. Summary and conclusions 

In this paper we present the data of a NIR photometric 
survey in the HDFN and Groth fields. Combining this data 
with the panchromatic sets available we have extracted a 
ii'-band selected sample up to K = 18.5 in a combined 
area of ~ 0.27 deg^. Additionally, we made use of the 
deep NIR public survey conducted by ESQ in the CDFS 
(172 arcmin^) to complement our relatively shallow sur- 
vey and increase the surveyed area. We have derived high 
quality photometric redshifts for the whole sample, taking 
advantage of the excellent coverage of the SED. X-band 
galaxy NCs have been derived in the three fields covering 
the range 16 < fsT < 18.5. The comparison to a compilation 
of shallow wide-area and deep pencil-beam surveys shows a 
general good agreement in spite of the considerable scatter, 
a consequence of the cosmic variance. We have studied the 
redshift distribution of the NCs finding that the 90% of the 
galaxy counts up to i^T ^ 18.5 come from the low-redshift 
population (z < 1). Additionally, we have measured the ef- 
fects of cosmic variance by comparing the redshift binned 
NCs in our three fields, finding that field-to-field differences 
reach 40% at a certain magnitude, with consequent impact 
on the NC distribution. 

We have provided evidence of an average decreasing 
trend in the slope of the if-band NCs, ~ 50% in the 15.5 < 
K < 18.5 range, evolving from dlog N/ dm = 0.6 — 0.3. The 
comparison to 100 synthetic catalogs from the Millennium 
simulation reveals that cosmic variance in areas ~ 0.25deg^ 
leads to significant scatter in the observed fiattening rate 
of the slope. 

We have studied the composition of the NCs in terms 
of LF building blocks, concluding that the change in the 
slope with observed magnitude can be summarized in three 
main regimes. Firstly, at bright magnitudes, the classical 
Euclidean regime [dlog N/ dm — 0.6) is dominated by low 
redshift M* galaxies. Then, at intermediate magnitudes, 
the transition regime is dominated by the LF at the redshift 
that maximizes 4'*^f- Here the slope decreases rapidly 
around the apparent magnitude of M* at that redshift. 
Finally, at faint magnitudes, the "a regime" is populated 
by galaxies at the faint end of a combination of LFs. Here 
the slope asymptotically approaches a minimum value at 
^-0.4(l-|-a). The value of a will typically be given by the 
LF at the maximum of 
larger faint end slope. 

We have explored the evolution of the fsT-band NCs by 
deriving LFs in the observed if-band as a function of red- 
shift (at a mean rest-frame wavelength 2.l6fim/{l + z)), 
and complementing our data with optical and NIR LFs 
from the literature. In terms of the multi-wavelength LFs, 
we find that the fiattening of the slope is the consequence of 
a prominent decrease of the characteristic density 4)*j^ ^j^^ (~ 
60% from z = 0.5 to z = 1.5) and the almost fiat evolution 
of M^^„bs (1^ compatible with M^^^^ = -22.89 ± 0.25). 
The combination of both effects forces a transition to the a 
regime at K ^ 17.5 that otherwise would have taken place 
1 — 2 magnitudes later. 

Our predictions for K >20, based on LFs from the lit- 
erature, suggest that the slope will continue to decrease 
steadily to dlog N/ dm <0.2 dominated by intrinsically 
faint galaxies from a mixture of LFs at z=1.5-2 and a mi- 
nor contribution from bright galaxies at z>2. However, if 
optical LFs at z > 3 present an even larger faint end, as 



or a close LF with a much 
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suggested by some authors, there might be a slope increase 
to a higher asymptotic hmit around K > 22.5. 
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